Strongly Interacting Luttinger Liquid and Superconductivity in an Exactly Solvable 

Model 
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A new family of exactly solvable one dimensional models with a hard-core repulsive potential 
is solved by the Bethe Ansatz for an arbitrary hard-core radius. The exact ground state phase 
diagrams in a plane 'electron density - on-site interaction' have been studied for several values of a 
hard-core radius. It is shown that superconducting phase and strongly interacting Luttinger liquid 
state are coexisted at a high electron density and unusually high value of repulsive on-site Coulomb 
interaction. 



The models of strongly correlated electrons with 
a bond-charge interaction which conserves a num- 
ber of double occupied sites, are simple examples 
of strongly correlated electron systems that exhibit 
superconductivityEl,El,B. The merit of these models is 
their complete integrability. The phase diagrams on a 
plane 'electron density - on-site interaction' have four 
phases, two of them exhibit off-diagonal long-range or- 
der (ODLRO) and thus are superconducting. The super- 
conducting phase is realized if the value of the on-site 
interaction less then the critical one U c and U c could be 
even positive at repulsive on-site interaction U. Other 
phases in which only single occupied and empty sites are 
presented called as U — * oo Hubbard state. In the case of 
a hard-core repulsive interaction between electrons with a 
hard core radius which exceeds a half of a lattice spacing 
the Luttinger liquid state transforms to strongly inte 
acting Luttinger liquid in a high electron density regior 
Strongly interacting Luttinger liquid is characterized by 
a large value of the critical exponent O for the momen- 
tum distribution function. At > 1 the residual Fermi 
surface disappears. At a high electron density when a 
hard-core repulsion interaction dominates, superconduct- 
ing phase and strongly interacting Luttinger liquid state 
do not coexist in the framework of the generalized t-J 
and Lai-Sutherland modelsa. The question arises: could 
the superconducting phase and strongly interacting Lut- 
tinger liquid state coexist simultaneously at positive and 
finite on-site Coulomb interaction? It turns out that an 
existence of the Fermi surface is not necessary for the 
superconducting phase. 

We shall consider a new family of integrable models 
and show that superconducting phase and strongly inter- 
acting Luttinger liquid states can coexist at a high elec- 
tron density or at small doping. This coexistence takes 
place at a repulsive on-site Coulomb interaction the value 
of which is large than a band width and depends on the 
value of the hard-core radius. In the modelsEl,0,B the hop- 
pings of single electrons on occupied states are forbidden, 
whereas the energy of electron pair is finite. In our mod- 
els we shall use the same hierarchy for the parameters 
of the interactions, the constants of interactions between 
single electrons are infinite and define a hard-core radius, 



the energy of electron pair is finite. We shall consider a 
new modification of a generalized one-dimensional Lai- 
Sutherland model for a study of a competition between 
strongly interacting Luttinger liquid state and supercon- 
ducting phase. The model Hamiltonian contains kinetic 
and interaction terms that combine those of the Hub- 
bard model and the Lai-Sutherland model. The model 
Hamiltonian includes two terms H. = Jihop + 'Hint 
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where Cj a and cj a are the creation and annihilation op- 
erators of fermions with spin er, a G {|, |}, L is the total 
number of lattice sites, < i,j > stands for neighboring 
sites, the projector Vi forbids two single electrons at dis- 
tances less than or equal to I ( I is measured in units of 
the lattice spacing parameter), t is the hopping integral, 
J is the constant of the exchange interaction. It is im- 
portant the Vi operator does not forbid doubly occupied 
lattice sites, as it takes place in the so-called U — ► oo 
Hubbard model or the t-J model. The last term in (2) is 
traditionally the most important term for the Hubbard 
model, the on-site Coulomb repulsion U separates the en- 
ergies of single and paired electrons states. The Hamilto- 
nian H conserves not only the total number of electrons 
N and also the number of single electrons with spin a 
Ni CT = y~)j*—i n ja{^ — n j-a) an d t ne number of localized 

electron pairs N 2 = J2j=i n j\ n 3l' N = Ha + 2N 2- 
In the case fcj:0 and J = the Hamiltonian (1), (2) 
is reduced toB and for fc=0, U — oo and J = t to 
the Lai-Sutherland modeH. For l> the Vi opera- 
tor is equivalent to additional two particle interactions 
between single electrons Yl r =i 52f=i ^r^ij^y+r; here 



according to (2) Ui+i = 



r ) with infinite U r parameters, 
-J. Using this representation 
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we can conclude that the kinetic term of the Hamilto- 
nian (1) is a particle-hole invariant: indeed applying this 
transformation c^ CT =>■ Cj a , Cj a ct ff , to the Hamiltonian 
(1), (2) we obtain J, (7) => H(t, J, U) + [7(L - N) 
. Due to a particle-hole symmetry the phase diagram is 
symmetrical with respect to a half filling. 

We examine the exact ground state phase diagram for 
the antiferromagnetic coupling J = t (we chose the hop- 
ping integral equal to unit then the coupling constants 
are dimensionless) and different values of the hard-core 
radius. The results of calculations are compared with 
the ones for J = - the simplest version of the model . 
Direct calculations show that the model (1), (2) is an ex- 
actly solvable one by the Bethe ansatz method and the 
set of the quasimomenta {kj}(j = 1,2, ...,Ni) satisfies 
the Bethe equations^ 
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and 



where P = % ^ s * ne momentum, Xj 

X a ( a — 1)2, ...,M) are the 'charge' and 'spin' rapidities 
M is the number of down spin single electrons. 

The eigenvalues and the magnetization are given by 



E = -2^cosfej + C/N 2 



3=1 



= i^N lCT -M. 



(4) 



(5) 



Let us introduce the partial electron densities: 
ni= Ni/L is the density of single carriers (Ni = 
53(7=1 i Nio-), n 2— N2/L is the density of localized elec- 
tron pairs. Clearly n = 111 + 2ii2, here n = N/L is the 
total density of electrons. 

Since the Bethe equations (3) we can calculate exactly 
the ground state phase diagram as a function of the elec- 
tron density and an on-site interaction for an arbitrary 
value of the hard-core radius. The density of localized 
pairs n 2 (or ni) can be calculated by minimizing the 
ground state energy per site £ = E/L for a fixed total 
density of electrons 



rQ 

£ = 2m -271-/ dAa(A)p(A) + Un 2 , 
J-Q 



(6) 



where a (A) = 



In the thermodynamic limit the Bethe equations re- 
duce to an integral equation of the Fredholm type for 
the function of the distribution of A, on the real axis 



p(A)+ / dA'K(A- A')p(A') = (1 -/m)a(A), (7) 



with the kernel being K(A) = ggzjg exp^A). 
The A— Fermi level denoted as Q controls the band fill- 
ing, the density of single electrons is defined by 



ni = 



Q 



dAp(A). 



(8) 



Q = corresponds to an empty subband of single carri- 
ers, ni = no for Q — > 00, here no = i+ gy 2 is a 'half-filled' 
density. Equations (6)-(8) are a consequence of the real 
solutions for 'charge' and 'spin' rapidities that describe 
the ground state of the system in the absence of an ex- 
ternal magnetic field. We calculate the critical exponent 
for a definition of the realization criterion of strongly 
interacting Luttinger liquid. We remind the momentum 
distribution function close to the Fermi momentum Uf is 
determined by the critical exponent 



(n fc ) 



K n k Fl 
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k F ), (9) 



where = — (1 



"and a = 2£ 2 (Q) is defined by the 



dressed charge C(A) according to the following integral 
equation C(A) + J Q Q dA'K(A — A')C(A') = 1-Zm. In the 
high electron density n > 2n2+n c (where n c is solution 
of equation 0(n c ) = 1 ) when the hard-core repulsive 
potential dominates the behavior of fermions is described 
as strongly interacting Luttinger liquid with > Id. 

We have focused on the calculation of the exact ground 
state phase diagram in the n— U plane for different val- 
ues of the hard core radius or I. First we consider pe- 
culiarities of behavior of the system using a simple ver- 
sion of the Hamiltonian (1),(2) when J = and then its 
transformation for J = 1. Due to a particle- hole sym- 
metry it is sufficient to discuss the phase diagram for 
n< 1. For J = the density of the ground-state energy 
(6) can be defined analytically £ = - 2i^ s i n ( t ^l_) + 

hU{n — 111), therefore a curve that separates a mixed 
region is defined according to the following equation 



■ cos 
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t/(m) = 4I sin (V^, — 

v Lf 7r VI — Zni/ 1 — Ini 

from -4 at m = to 4(1 + I) at n 
at ni = to 4(1 + I) at ni = n max = 1/(1 
maximal value U c = 4(1 + 1). The value of n c is equal to 

n c = (1 ' 
for I = 2, n 



, - — , . U variates 
U variates from —4 
I), hence a 



y/6~4:V2)/l (n c = 0.414 for I = l,n c = 0.207 
0.138 for 1 = 3). For I = 1 the complete 
phase diagram is shown in Fig. |l|. The lower region ( for 
U < —4 ) is characterized by only double occupied (solid 
circles) and empty (empty circles) sites so ni = and 
n 2 =n/2. For —4 < U < U(ni) we have a mixed region, 
the ground state includes both finite densities of single 
electrons (spheres with dot center) and localized electron 
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FIG. 1: Ground state phase diagram in the case J = for 
1 = 1. Dashed line corresponds to the modelcl (or / = 0), 
dotted lines separate metallic and insulator phases. The area 
of a strongly interacting Luttinger liquid state is denoted as 
SILL. 



pairs. These phases are superconducting, since the two 
particle correlation function < T]\rij > (here m — cLct^ 

) exhibits ODLRO0,! i.e. < r}\r}j >^ for | i — j |— > oo. 
At n>n c a curve f7(n c ) separates this phase on Luttinger 
liquid state at —4 < U < U(n c ) and strongly interact- 
ing Luttinger liquid state at U(n c ) < U < [/(ni).The 
strongly interacting Luttinger liquid state is denoted as 
SILL in the figures. Note, that strongly interacting Lut- 
tinger liquid state is realized at largest values of a re- 
pulsive on-site Coulomb interaction and a high electron 
density. Comparing the calculations for different I we 
can conclude that a hard-core repulsive interaction in- 
creases a region of the coexistence of strongly interact- 
ing Luttinger liquid and superconducting phase due to 
both a larger U c and smaller n c . For U > U(ni) and 
n < n max the ground state coexists of singly occupied 
and empty sites; dotted lines separate a metallic phase 
(at n < n max ) and an insulator phase (at n > n max ) with 
a gap Ae = U — U c . 

An exact solution of the problem enables to study the 
role of the exchange interaction on the behavior of a 
strongly interacted electron system. Let us consider a 
transformation of the exact ground-state phase diagram 
for J — 1, we restrict our consideration the case n < n . 

increases monotonically from | to ft — (3+2;)^ 2 

with the ni density. We should solve equation 8(n c ) = 1 
numerically calculating the dressed charge as a function 
of the electron density n! for arbitrary I; for example 
n c = 0.348 for I = 1, n c = 0.192 for I = 2, n c = 0.131 



FIG. 2: Ground state phase diagram in the case J = 1 and 
1 = 2- similar to that for Fig. |l|. 



for I = 3, n c = 0.1 for 1 = 4. According to the numerical 
results obtained the critical density n c is less than no. For 
/ = 2 the ground state phase diagram is given in Fig. g. 
All electron states: empty, single occupied and doubly 
occupied sites are presented simultaneously in a mixed 
region (a closed region in figure 2). For n c < n < no two 
branches of curves separate Luttinger liquid state and 
strongly interacting Luttinger liquid, that is realized be- 
tween these branches. Comparing the phase diagrams for 
J = 1 (Fig. g) and J = (Fig. 0) calculated for the same 
value of I we can conclude that the exchange interaction 
decreases the region of the coexistence of the supercon- 
ducting phase and strongly interacting Luttinger liquid 
state (n c and U c decrease slightly) . U c increases with an 
increasing of hard-core radius. 

In summary, we have presented a soluble generalization 
of the Lai-Sutherland model, having the nontrivial Lut- 
tinger liquid behavior. The exact solution was obtained 
by means of the nested Bethe Ansatz. We have derived 
the exact ground-state phase diagram; the latter exhibits 
an unusual phase state in which superconducting phase 
and strongly interacting Luttinger liquid state coexist. 
This phase is realized at high electron density and posi- 
tive values of the on-site Coulomb interaction. The maxi- 
mum critical value U c realized in the mocLeLis higher than 
that of all other exactly solvable modclsojj. This is im- 
portant because higher values of U c expands the region of 
coexistence of superconducting phase and strongly inter- 
acting Luttinger liquid state. It has been shown that the 
presence of the Fermi level is not necessary for realization 
of superconducting phase. The results of calculations of 
one dimensional models do not allow direct application 
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to the real 2D and 3D systems. Nevertheless one can 
assume that real high-T c superconductors belong to the 
family of strongly interacting Luttinger liquid described 
above. 
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